Thermodynamic phase transition of UH 3 : Role of electronic strong correlation 
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The electronic structure and thermodynamical properties of uranium trihydrides (a-UH3 and /3- 
UH3) have been studied using first-principles density functional theory. We find that inclusion of 
strong electronic correlation is crucial in successfully depicting the electronic structure and thermo- 
dynamic phase stability of uranium hydrides. After turning on the Hubbard parameter, the uranium 
5/ states are divided into well-resolved multiplets and their metallicity is weakened by downward 
shift in energy, which prominently changes the hydrogen bond and its vibration frequencies in the 
system. Without Coulomb repulsion, the experimentally observed a— 1/3 phase transition cannot 
be reproduced, whereas, by inclusion of the on-site correlation, we successfully predict a transition 
temperature value of about 332 K, which is close to the experimental result. 

PACS numbers: 71.27.+a, 71. 15. Mb, 63.20.D- 



Since the report of UH3 as one ferromagnetic (FM) 
uranium compounds^— , it has attracted lots of attentions 
in nuclear fuel field in the past decades. In solid state, 
uranium hydride exists mainly in the form of cubic trihy- 
dride (Pm3n; No. 223)^ with two allotropes, i.e., a-UH 3 
and /3-UH 3 . a-UH 3 [Fig. 1(a)] is an intermediate state 
from orthorhombic a-U saturated with H to ^-UH^— . 
It is the metastable low temperature phase with two UH3 
formula units in one unit cell. The uranium and hydrogen 
atoms in this structure occupy the 2(a) (0, 0, 0) and 6(c) 
(1/4, 0, 1/2) sites, respectively. /3-UH 3 [Fig. 1(b)] is the 
stable high temperature phase with eight UH 3 formula 
units in one unit cell. The uranium atoms in this struc- 
ture occupy the 2(a) (0, 0, 0) and the 6(c) (1/4, 0, 1/2) 
sites, and the hydrogen atoms occupy the 24(k) (0, 0.156, 
0.313) sites. Mulford et al& found through X-ray diffrac- 
tion measurement that the transition temperature of a- 
UH 3 to /3-UH 3 is between 373 K and 523 K, while Genos- 
sar et alJ- reported that a-UH3 transits irreversibly into 
/3-UH 3 above room temperature. At even higher temper- 
atures (above 673 K), uranium hydride reversibly remove 
the hydrogen. This property makes uranium hydrides 
convenient starting materials to create chemically reac- 
tive uranium powder along with various uranium carbide, 
nitride, and halide compounds. In contrast to the above 
mentioned experimental measurements on temperature- 
induced a— >/3 phase transition of uranium trihydride, to 
date the theoretical exploration of this prominent phase 
transition is still totally lacking in the literature. Con- 
sidering the combined fact that (i) UH 3 denotes a typical 
prototype among various kinds of actinide hydrides and 
(ii) the actinide hydrides play extremely important role 
in nuclear fuel design (as well as stockpile) in which the 
thermodynamic stability is a substantial factor, there- 
fore, a revealing theoretical study on the ground-state 
properties and a— >/? phase transition of UH 3 from first- 
principles quantum mechanics is highly needed for the 
relevant important industrial applications. 

From basic point of view, it can be visualized that 
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FIG. 1: (Color online) Crystal structures of (a) a-UH3 and (b) 
/?-UH3, where the large and small balls denote the uranium 
and hydrogen atoms, respectively. 



many physical and chemical properties of UH3 are closely 
related to the quantum process of localization and der- 
ealization for partially filled uranium 5/ electrons. Mod- 
eling of the electron localization/delocalization, and thus 
any hydrogenation process involving uranium, is a com- 
plex task. Conventional density functional theory (DFT) 
schemes that apply the local density approximation 
(LDA) or the generalized gradient approximation (GGA) 
underestimate the strong on-site Coulomb repulsion of 
the uranium 5/ electrons and consequently fail to cap- 
ture the correlation-driven localization. Therefore, the 
uranium 5/ electrons require special attention in trying 
to gain any microscopic insight into the thermodynami- 
cal stability of UH3. Up to now several approaches, the 
LDA/GGA+t/, the hybrid density functional of (Heyd, 
Scuseria, and Enzerhof) HSE, the self-interaction cor- 
rected local spin-density (SIC-LSD), and the Dynami- 
cal Mean-Field Theory (DMFT) , have been developed to 
correct the pure LDA/GGA failures in calculations of ac- 
tinide compounds. Among these techniques, the effective 
modification of pure DFT by LDA/GGA+J7 formalisms 
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has been confirmed widely in study of uranium (as well 
as its neighbor, plutonium) oxides^—. By tuning the ef- 
fective Hubbard parameter in a reasonable range, the an- 
tiferromagnetic (AFM) Mott insulator features of these 
oxides were correctly calculated and the atomic struc- 
tural parameters as well as the electronic properties are 
well in accord with experiments. Inspired by this series 
of successful calculations on AFM uranium oxide, as well 
as motivated by the fact that uranium hydride is also a 
magnetic ordered (FM instead of AFM) compound and 
thus a static mean-field treatment like LDA/GGA+f/ is 
expected to reliably catch its some key features in both 
a and /3 phases, in the present work, we investigate the 
ground-state electronic properties and lattice dynamics 
of the two allotropes of UH3 using the LDA+ U formal- 
ism. For comparison the pure LDA calculation is also 
performed. The most insightful result of our investiga- 
tion is that it is essential to take the 5/ on-site electronic 
correlation into account for theoretically reproducing the 
experimentally observed a— S-/3 phase transition of UH 3 at 
finite temperature. Also, the prominent changes in the 
hydrogen bond and the consequent hydrogen vibration 
(reflected by the optical branches in phonon spectrum) 
in UH3 by the inclusion of on-site Coulomb repulsion of 
the uranium 5/ electrons is for the first time highlighted 
in this paper. 

The calculations are performed using the projector- 
augmented wave (PAW) method of Blochl 11 , as im- 
plemented in the Vienna ab initio simulation program 
(VASP) 12 . For the plane-wave set, a cut-off energy 
of 500 eV is used. The hydrogen Is and uranium 
6s 2 6p 6 6d 2 5/ 2 7s 2 are treated as valence electrons. The 
strong on-site Coulomb repulsion amongst the localized 
uranium 5/ electrons are accounted for by using the for- 
malism formulated by Dudarev et al^. In this scheme 
the total energy functional is of the form 



Elda+u — Elua 



U -J 



£[Tr^-Tr(pV)], (1) 



where p a is the density matrix of / states, and U and 
J are the spherically averaged screened Coulomb energy 
and the exchange energy, respectively. Here the Coulomb 
U is treated as a variable, while the exchange energy is 
set to be a constant J=0.51 eV. This value of J is in 
the ball park of the commonly accepted one for uranium. 
Since only the difference between U and J is significant^, 
thus we will henceforth label them as one single parame- 
ter, for simplicity labeled as f/ e ff, while keeping in mind 
that the non-zero J has been used during calculations. 
The Brillouin-zone (BZ) integrations are performed using 
9x9x9 and 7x7x7 Monkhorst-Pack 14 special fc-points 
for cv-UHj and /3-UH3, respectively. The geometries are 
optimized until the forces are less than 0.02 eV/A, and 
the total energy is relaxed until the difference value is 
smaller than 10 -5 eV. 

Through calculating the total energy dependences on 
U e ff of a-UH 3 and /?-UH 3 in nonmagnetic, FM, and AFM 
phases, we find that the FM phase is the most favorable 
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FIG. 2: (Color online) Spin-resolved total DOS for the two al- 
lotropes of UH3 calculated within LDA and LDA+ U ( U e g=4 
eV) formalisms. Partial DOS of U 5/ and H Is orbitals are 
also shown. The fermi level is set to be zero. 



state with U e g from to 7 eV, which is well consistent 
with the experiment. The dependence of the lattice pa- 
rameters on f/ c ff for FM phase of CV-UH3 and /3-UH3 
demonstrates that in a phase, the lattice parameter is 
close to the experimental value when U e g is in range of 
4-6 eV, while in /3 phase, the range is in 4-5 eV. Consid- 
ering the two allotropes, therefore, we choose the value 
of U c s to be 4 eV in our following study for the two 
phases. Our calculated lattice constants for a (/?) phase 
of UH 3 at f/cff =0 and 4 eV are 4.001 (6.382) and 4.128 
(6.622) A, respectively. Compared with pure LDA, our 
LDA+f/ gives closer values with respect to the experi- 
mental values of 4.160^ and 6.639^ A for a and j3 phase, 
respectively. 

Our calculated spin-resolved total density of states 
(DOS) and the partial DOS for uranium 5/ and hydro- 
gen Is states of the two allotropes are presented in Fig. 
[U The upper panels show the pure LDA results, while 
the lower ones give the LDA+ U results. From Fig. [5] we 
can observe evident effects of strong correlation on the 
/ states. For both two phases of UH 3 , after turning on 
the Hubbard U parameter, clearly, the uranium 5/ states 
are divided into more well-resolved peaks compared with 
pure LDA, and the energy distributions of these multi- 
ple 5/ peaks are much localized and separated. As a 
result of this multiple peak splitting, for a phase, the 
uranium 5/ states in the range of —1.0 to 2.0 eV within 
LDA are extended to lie in the range of —2.5 to 5.0 eV 
in the LDA+ U formalism. For (3 phase, similar effects 
of strong correlation can be found. Moreover, the occu- 
pation of both spin-up and spin-down electrons at the 
Fermi level are lowered after taking into account the on- 
site Coulomb repulsion. However, no insulating band gap 
are opened by the inclusion of Hubbard parameter and 
UH3 still remains its metallic nature as observed in the 
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TABLE I: Bader effective atomic charges and volumes of a- 
UH 3 and /3-UH 3 in the LDA and LDA+ U (U cB = 4 eV) for- 
malisms. 

Allotrope Methods Q(U : ) Q(U n ) Q(H) V(UQ V(U n ) V(H) 
a-UH 3 LDA 12.527 1.490 17.913 4.703 

LDA+ U 12.425 1.525 18.847 5.415 

/3-UH 3 LDA 12.503 12.434 1.503 18.290 17.465 4.904 

LDA+U 12.452 12.370 1.535 18.774 19.283 5.706 



Due to the decrease in metallicity of uranium 5/ states 
and their energy downward shift towards the hydrogen Is 
orbital by taking account of the electronic strong on-site 
correlation, the ionicity of the hydrogen bond is promi- 
nently enhanced. To analyze the ionicity of the two al- 
lotropes of UH3, results from the Bader analysis^! are 
shown in Table HI The charge ( Q ) enclosed within the 
Bader volume ( V) is a good approximation to the total 
electronic charge of an atom. Note that although we have 
included the core charge in charge density calculations, 
since we do not expect variations as far as the trends 
are concerned, only the valence charge is listed. From 
the results we calculated, it is seen that the ionicity of 
the two allotropes within LDA+ U are more evident than 
that within pure LDA. For a phase, the electrons trans- 
ferred from each U atom to H atoms are 1.473 and 1.575 
within LDA and LDA+ U formalisms, respectively. For 
P phase, electrons transferred from two inequipotential 
U atoms (Ui and Un) to H atoms are 1.548 and 1.630, 
respectively, in LDA+ U formalism. Whereas, in LDA 
formalism the values are 1.497 and 1.566 respectively. 

To illustrate the main point of this paper, i.e, the effect 
of electronic strong on-site correlation on the thermody- 
namic properties of UH3, we have calculated the phonon 
dispersion and the Helmholtz free energy F of two al- 
lotropes with and without including Coulomb repulsion 
of 5/ electrons. Phonon frequency calculations were car- 
ried out using the Hellmann-Feynman theorem and the 
direct method^. For a and P phases, we adopted 2x2x2 
supercells containing 64 and 256 atoms with 3x3x3 and 
lxlxl Monkhorst-Pack fc-point meshes in the BZ inte- 
gration, respectively. The forces induced by small dis- 
placements were calculated within VASP. The amplitude 
of all the displacements is 0.03 A. 

The calculated phonon curves along some high- 
symmetry directions in the BZ together with the phonon 
DOS are plotted in Fig. 3. Clearly, the effect of elec- 
tronic strong correlation on phonon dispersion can be 
observed by comparing results from LDA and LDA+ U 
approaches. Compared with pure LDA calculation, the 
optical branches calculated within LDA+ U shift down 
by around 2.5 THz and 2.3 THz for a-UH 3 and /3-UH 3 , 
respectively. However, the acoustic branches have no ev- 
ident changes for both phases. Due to the fact that ura- 
nium atom is much heavier than hydrogen atom, then the 
optical branches denote the hydrogen vibration while the 
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FIG. 3: (Color online) Phonon dispersion curves (left panel) 
and corresponding phonon DOS (right panel) of (a) q-UH 3 
and (b) /3-UH 3 . The black curves are computed within the 
LDA while the red (grey) curves the LDA+ U formalism. 



acoustic branches come from the uranium vibration. As 
a result, a large gap between optical and acoustic modes 
can be observed. The key point revealed in Fig. 3 is 
that the lattice dynamics behaviors are prominently in- 
fluenced by electronic strong on-site correlation, which 
has changed the U-H bonding nature. This influence on 
lattice dynamics turns out, as analyzed below, to be fun- 
damental to reproduce the experimentally observed rela- 
tive thermodynamic stability sequence of a and P phases 
ofUH 3 . 

For the sake of determining phase transition tempera- 
ture of a-UH3 and /3-UH3 at ambient condition, we have 
calculated the Helmholtz free energy F in the LDA and 
LDA+ U formalisms. This quantity at volume V and 
temperature T can be expressed as 

F(V, T) = E(V) + F vih (V, T) + F clc (V, T), (2) 

where E(V) is the K band energy, F v ib(V, T) is the 
phonon vibrational free energy, and F e \ e is the thermal 
electronic contribution to the free energy. The phonon 
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FIG. 4: (Color online) Temperature dependences of the 
Helmholtz free energies within LDA (the black curves) and 
LDA+?7 (the red/grey curves) formalisms. The inset gives 
closer view of the crossing point within the LDA+ U. 

vibrational free energy in the quasiharmonic approx- 
imation can be calculated from phonon DOS g(oJ, V) 
as F vih (V, T)=k B T r Q °° g(u, V) lnpsinht^)]^. The 
thermal electronic contribution can be described 
as F C \ C (V, T)~E C \ C — TS e i e with electronic entropy 
S ele (V, T)=-k B J n{e, V)[fhx f+(l-f) m(l-/)]cfc, 
where / is Fermi-Dirac distribution and n is the elec- 
tronic DOS. The chemical potential at temperature T is 
fixed by conservation of total valence electron number 
during thermal excitation. It turns out that the value of 
the thermal electronic contribution is very small, only 
one-tenth of the phonon vibrational free energy. 

We have calculated and plotted the temperature de- 
pendences of the Helmhotz free energy for the two UH3 
phases within LDA and LDA+Z7 formalisms (Fig. 0]). A 
significant difference can be clearly observed between the 
two formalisms. Without including the on-site Coulomb 
interaction, the Helmhotz free energies of two UH3 phases 
have no crossing point, and the a phase is always more 
thermodynamically stable than the (3 phase in a wide 



temperature range of to 1000 K, as shown in Fig. 
2] This is contrary to the experimental observations. 
Whereas, after including the electronic strong correla- 
tion, the calculated Helmhotz free energy curves of the 
two allotropes cross at the temperature of 332 K, which 
is well consistent with the experimental value for a— >/3 
phase transition. Therefore, we arrive at that while the 
pure LDA calculation totally fails to reproduce the exper- 
imentally determined thermodynamic phase transition of 
UH 3 , the LDA+ U calculation can well depict it. This is a 
good news for the science of actinide hydrides. Needless 
to say, more advanced many-body techniques, such as 
DMFT, which can account for some spin fluctuations at 
finite temperatures, will further improve the calculations 
and we would like to leave it for future consideration. 

In conclusion, the ground-state properties of two al- 
lotropes of UH 3 have been comparatively studied within 
the LDA and LDA+ U formalisms. After switching on 
the on-site Coulomb interaction, the uranium 5/ states 
split into the multiple well-resolved peaks that shift up- 
ward or downward with respect to the Fermi energy, 
which reduces the metallicity of the 5/ states and en- 
hances the ionicity of the hydrogen bond in the system. 
As a result, the calculated optical branches in the phonon 
dispersion of two allotropes have been found to shift down 
by about 2.5 and 2.3 THz for a-UH 3 and /8-UH3, re- 
spectively, by turning on the Hubbard parameter. Our 
theoretical Helmhotz free energy curves have shown that 
while the pure LDA fails to describe the a— >(3 phase tran- 
sition, the LDA+ U calculation gives the transition tem- 
perature of 332 K, which well lies within the experimen- 
tally observed range. These results clearly indicate that 
the electronic strong on-site correlation plays an impor- 
tant role in the uranium hydride systems. 
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